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Abstract. Synthesizing algorithms of decentralized adaptive control for a class of interconnected systems 
with nonlinear local dynamics and uncertainty are considered, ensuring the stability of isolated 
subsystems of the entire system, taking into account interactions. When connected to a stable isolated 
subsystem, the system's stability may be violated and the objective function will cease to be a Lyapunov 
function for a system with connections. Very strong connections can significantly change the dynamics of 
the system, and the behavior of subsystems in an isolated state will be radically different from the 
behavior in a coupled system. Therefore, the law of control of a complex system, which ensures only the 
stability of isolated subsystems, cannot guarantee the stability of the entire system considering 
interactions. At the same time, to solve the problem of taking into account the interactions of subsystems, 
an adaptation algorithm based on the speed gradient method was introduced into the general control law. 
The results of the numerical analysis confirmed their effectiveness, which allows them to be used in 
solving applied problems of synthesizing systems for decentralized adaptive control of technological 
objects. 
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Introduction 

According to theoretical point of view, completely centralized systems in control theory are 
considered ideal, allowing for the greatest effect. However, the implementation of a single controller that 
measures the entire set of state variables and generates the entire set of controls in a complex system is 
fraught with great difficulties and requires high costs. In practice, a system is often built in such a way 
that, due to its design features, it is divided into localized subsystems that have a certain degree of 
autonomy in the sense that the subsystem can function and be controlled not only by a system-wide body 
but also independently by a local regulator. If the original system allows for division into a number of 
interacting systems, then the decentralized control algorithm for each subsystem generates a control action 
based on information about the state of only this subsystem, without taking into account the state of other 
subsystems [1-8]. In this case, control synthesis is usually preceded by a simplification of the problem, 
which consists in discarding connections between subsystems. Therefore, the problems of adaptive 
decentralized control of interconnected systems occupy an important place in the theory of control of 
complex systems. 

1. Problem formulation 

Let us consider a control object (CO), consisting of interacting subsystems, the dynamics of each 
N of which is described by the following equations: 

, = f(x) +B,(%,,u,) +h,0); 1=1,2,...,N, (1) 


; th 


where x, € R” — state vector of the i” subsystem, u, € R” — input vector of the i” subsystem, 


N N 
n= hs me me X =COI(X,,Xy5..5Xy) € R" — cumulative system state vector, f;(-), b,(-), h;(-) 
i=1 i=1 


— continuously differentiable vector functions. 
Thus, CO (1) is a complex system consisting of N interconnected nonlinear subsystems, where the 
f,©,. 6,©) functions characterize the own dynamics of thei ” subsystem, and the h,(-) functions 


describe the connections (interactions) between the subsystems. We will assume that all the uncertainty of 


the subsystems lies in the connection functions h,(-) satisfying the following inequality: 


N 
h; aye * 
j=l 


Setting the task: to build a decentralized control algorithm: 
uj =U; (x;,9), (3) 


.>0,i=1,2....,N. (2) 


i,j 


6. =9,(x,,6), (4) 
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where 6. € R’ — adjustable parameters ensuring the limited trajectories of the closed-loop system 
(1), (3), (4) and achieving the control goal: 

Q;(x;(t)) > 0 npu too, (5) 

where Q.(x;) 20 — given smooth objective functions, i =1,2,...,N. 

The task is significantly simplified if we discard connections between subsystems and consider the 
CO (1) as a set of isolated subsystems, Le. 

x; = f,(%,) +b,(%;,,u;), #=1,2,...,N. (6) 

Then control algorithms u; (x;) that ensure the fulfillment of quadratic type inequalities for local 
objective functions [2, 3, 9-12] 


a,|x[ $O,(%,) $ a,|x,) 


2 


2 


VO.(%)|< Bx. (7) 
O14}, aro) SPI 
where QO:(%), cy =VO,(x;)" (F,%;) +B: (%;.u; (4); ys A,B, P, >0, 1=1,2,...,N, 


guarantee exponential stability of the entire CO. 

However, the presence of interactions between subsystems qualitatively changes the behavior of a 
complex system (1), and the fulfillment of inequalities (7) for each isolated subsystem obviously does not 
mean the stability of the entire system (1). Indeed, for a system without interactions (6), the derivative of 
the objective function AA has the form: 

O; (x;54;) = O;(5)" (F%;) +B, (,4;)s 

and for a system with connections (1): 

O,(x;4;,) = VO)" (F.(%)) +B, ju) + h,00), 

i.e. the stability of each subsystem depends on other subsystems. Thus, when connected to a stable 
isolated subsystem of connections, stability may be violated and the objective function will cease to be a 
Lyapunov function for a subsystem with connections. 

Very strong connections can significantly change the dynamics of the system, and the behavior of 
subsystems in an isolated state will be radically different from the behavior in a connected system [3, 12- 
14]. Therefore, the control law for a complex system (1), which ensures only the stability of isolated 
subsystems, cannot guarantee the stability of the entire system taking into account interactions. 

2. Construction of a decentralized adaptive control law 

Let assume that each i” isolated subsystem (6) is exponentially stabilized, i.e., control laws u, (xy 
i=1,2,...,N are known that satisfy inequalities of type (7). 

Let’s assume that the functions U; (x;) can be represented as follows: 


ui; (X;) =U; (%j,8;), 1=1,2,..4N,, 
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where @, —some constant parameters (for example, gains involved in the control law). 

Then the control algorithm for the system with interactions (1) is chosen as following: 

u, =U;(%;,6,), i=1,2,..,N, (8) 

Where the function U;(-) — the same as before, however, instead of the known constant parameters 
@, in the control law (8), their settings 6. are used, calculated by the adaptation algorithm, based on the 
speed gradient method. 

To synthesize the adaptation algorithm, following the speed gradient method, the total derivative of 
the objective function Q,(x;) along the trajectories of the system was calculated (6), (8): 

O,(x;,0;) = VO)" (Fi) +B (%4.Ui(%,8,))) 

Then the velocity gradient looks like this: 

V 4 Qi(%1,8,) = Vg (O:(4)" BU (4,8), 

and the adaptation algorithm looks like: 

6 =e (VO. (x,)" b,(x,,U;(%;,4,))), P= Ll: (9) 

Taking into account the constructed adaptive control law (8), (9), we will consider new Lyapunov 
functions V,, (x,,4,) for subsystems (6), depending on the adjustable parameters 6. eR’ and having the 
following form: 

Vi (ir 4) = O10) + 5-0-0)", -8). (9) 


3. Properties of a synthesized control system 
Let the following conditions be satisfied for CO (1): 


a) for each i" isolated subsystem (6) there are smooth functions U; (x,,0;) , vectors and positive 


numbers @,,,@,;,,,,,1 =1,2,...,N such that the following inequalities are satisfied: 


Oy ral < O(x%,) Sa, |x, ; 
VO(x)|< Bx, 

OAR S,)) «ea; = —O.O): 

where O,(%,,9,)|,,-0.¢+,,0) = O:(%-8) = VO,(%,)" (F,(4;) + B(%,,U; (%,.8,))), — condition of 


stabilization of the i” isolated subsystem; 


° 


(10) 


° 


b) the complete derivatives of the objective functions Q;(x;) are convex in 6, ie. for any 6.61, Xj, 
the convexity conditions are satisfied—the following inequalities: 
O;(x;,91) = O;(x;,8,) 2 (,-9, ) Vg (O;(%.9,)), (11) 


c) interaction functions h,(-), i =1,2,...N satisfy condition (2) and the additional condition: 
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N 
B; +> Sy < 22,0; /VN, i=1,2.,...,N. (12) 
jl 


Then the decentralized adaptive control law: 
=U(x,,0), 

“i L (x; L ) (13) 

0, — —7ViV 5 (VO'b, (x,,U;(%,.8,))), 

where 7, > 0,1 =1,2,...N, ensures the limited trajectories of the closed-loop system (1), (13) and 


the achievement of the control goal (5). 
Let us define the Lyapunov function for a complex system (1) as the sum of the Lyapunov functions 


of subsystems (9): 
V, (x, 0) = ay (x,,0) = dia (x,) + 34- 0)" (8, -8)} 
and calculate its total derivative. By virtue i _— (1), (13), we write: 


Vv, = Sf va! (fea +b th 6o.8)-+h)+2O -4y'6} - 


N aw aw 

= D1 VO? (F.%,)+8,,U, 8) + VOTH) - 6, -8)'V,(VO7H)}. 
i=l : 

Using the convexity condition, we have: 


V, = > {6 O,(%,,8) - 6 -8)'V,(0,) + VOlh(o} < <>{0,0%,,.0)+ VO}. 


i=l 
Then from the condition of stabilizability of isolated subsystems and restrictions on the functions of 
relationships h,(-) (2) and target functions Q,(-) (10) it follows that: 


N N N N 
> —> p,2,(x,) + DB [,| bar |x, | = —>) ap; [x,| ‘ 
i=l i=l j=l i=l 


Applying to the last expression the well-known scalar inequality 2ab<a*/c+b’c valid for c>0, 


7 


we have: 


ij 


DD} Ait . file cae 


x, AL re} 


N N, 
- rE prs tad feed 
i J i i 
Let us introduce the following notation: 


aXe Da}-e, 


By choosing h,(-), satisfying inequalities (2) and (12), we obtain that 7; > 0. As a result we have: 
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; N 5 N n; 
V,5-) mia] =—>,—- O@). (14) 
i=l i=l 22; 
where 77,/@,; >0. 
Thus, V, (x,0) <0, whence it follows that Vy, (x(t),0(t)) < V, (x(0),0(0)), ice. V, (x(t), A(t) , does not 


increase and is bounded from below, and therefore has a finite limit at tf oo. 
Integrating (14) over the interval [0, t], we obtain the following inequality: 


0< YAL,()) <V, 0,40) <V, 20), 40) - )™ [O,4%,(s)0ds, 


i 2i 0 


whence it follows that } Q.(x,(t))dt i=1,2,...,N , exist and are finite. 
0 


The boundedness of the trajectories x(t), O(t) of the system (1), (13) follows from the growth 
condition (10) of the functions Q,(-),i=1,..., N. 

Further, taking into account the boundedness of the right-hand sides of the closed-loop system (1), 
(13), it can be shown that all x,;(¢) are uniformly continuous, and, therefore, all functions Q;(x;(t)) are 
uniformly continuous and integrable on [0,0) . Then, using Barbalata’s lemma [3,14—17], we obtain that 
Q;(x;(t)) >0 for toa. 


The adaptation algorithm in (13) can be obtained using the so-called diagonal dominance principle 
[3, 18]. Note also that condition (a) ensures the achievement of the control goal (CG) in each i” isolated 
subsystem, and condition (c) determines the permissible degree of interconnection of subsystems. 
Connections that satisfy condition (b) can be called connections of a homogeneous type. 

Consideration of interactions of a heterogeneous type, subject to weakened inequalities: 


N 
h(x) as yc; 
j=l 


x;|-+d;, (15) 


where d; >0,i=1,...,N, allows you to take into account the influence of limited disturbances. In 


this case, we cannot guarantee that the original control goal (4) will be achieved. However, as will be 
shown below, it turns out to be possible to provide a weakened control center that ensures the 
convergence of system trajectories into a certain limited limit set in the state space of the system. It should 
be noted that the adaptive control law (13) cannot ensure the achievement of even the specified weakened 
goal [17-19], therefore the control algorithm (13) must be coarsened in various ways: for example, by 
introducing a dead zone or, as shown in the following theorem, by introducing negative feedback 
adaptation algorithm. 
Let the interaction functions satisfy condition (15), and the additional condition: 


N 
BP+ dE} <2p,a,/VN, i=1,2,..,N. (16) 
j=l 


Then the decentralized adaptive control law: 
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U; =U,(x;,6,), 
6, =-AAV (VO! (xU,(x;,6)))- 0,8}, 


where y; >0, a; > 0 — constant gain factors, i=1,2,...,N, ensures limited trajectories of the closed- 


(17) 


loop system (1), (17) and achieves a weakened control goal: 
lim Q,(x;(1) < A;, (18) 
t>o0 


where A, = 0, i =1,2.....N. 


Let us define the cae function for system (1), (15) as the sum of Lyapunov functions (9): 


V5 (9)= 2 (30) = Y1OK)+ 5 8- 0,)' (0, -6,)}. 


1 


and calculate its total derivative due to system (1), (17): 


N 

V, = (VO! (F.(x;) + b,(%;,U (3%, 8) + VOF A, (x;) - (8, - 8)" Vg (VO; b;) - (6,-06,)' 0,6}. 
i=1 

Then, agreeing with the above method, we will have the following: 


V, = ¥105.8)- (6, -8,)'V,(Q)+VOrh (x) - 
~(6-0)" 08} < 10,4.) +VO"h(x)- 6-0)" oA). 


Then from the condition of stabilizability of isolated subsystems and restrictions on the 
interconnection functions h,(-) (15) and objective functions Q,(-) (10) it follows that: 


< —D PQ, x|(; 


N 
= -) QP; 
i=l 


x,|+d,)-@ -~0)"oO}< 


y 


N 
ld; ic, @, -6,)' 6. 
i=l 


272 
IIpu o, (6, =6," 425 6;| Pw Bedi|a|< la) fie 
i 
It can be written as follows: 
eB bel . - 
i=l 2 2 
N 2 
¥] 2 o)-Z 0 -af +A. 2p 
Denoting p; = min{ ! , 0;y;}. Then: 
Ay; 
N Bd? at. 
¥,33- “PQ +5 - 0)" (6,- ay a Sar 
i=] i 
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Introduce denotation 
A 

A; = p;, B; 

Then 


N N 
V, = >.V, = D{-AV, +B}. 
i=l i=l 
Integrating the resulting differential inequality, we obtain: 


O,(x;(t)) < V,(x; (1)0,(t) < V,(x,(0)0,(0)) yy" i i=1,2,...,N, 


whence it follows that the trajectories x,(f) , 6, (t) are bounded and, in addition, 
lim Q;(x;(t)) < A;, 
t>o0 


where A, ae L= 1, 2eeug lV « 
B; 
Conclusion 
Obtained results make it possible to ensure not only the stability of isolated subsystems, but also the 
stability of the entire system, taking into account interactions. At the same time, to solve the problem of 
taking into account the interactions of subsystems, an adaptation algorithm based on the speed gradient 
method is introduced into the general control law. 
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